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Abstract: With the Deep Underground Neutrino Experiment (DUNE) as an example, we 
show that the presence of even one sterile neutrino of mass ~1 eV can significantly impact 
the measurements of CP violation in long baseline experiments. Using a probability level 
analysis and neutrino-antineutrino asymmetry calculations, we discuss the large magnitude 
of these effects, and show how they translate into significant event rate deviations at DUNE. 
Our results demonstrate that measurements which, when interpreted in the context of the 
standard three family paradigm, indicate CP conservation at long baselines, may, in fact 
hide large CP violation if there is a sterile state. Similarly, any data indicating the violation 
of CP cannot be properly interpreted within the standard paradigm unless the presence 
of sterile states of mass 0(1 eV) can be conclusively ruled out. Our work underscores 
the need for a parallel and linked short baseline oscillation program and a highly capable 
near detector for DUNE, in order that its highly anticipated results on CP violation in the 
lepton sector may be correctly interpreted. 
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1 Introduction 

A major goal of present and future long-baseline neutrino oscillation experiments is to 
establish that leptons violate CP, or else to place a stringent upper limit on any such 
violation. (For recent status reviews see [1-12].) Our thinking about these experiments 
usually assumes the standard neutrino paradigm, in which there are just three neutrino 
mass eigenstates separated by just two independent mass-squared splittings, three mixing 
angles Oij, and just one CP-violating phase Jcp relevant to oscillation. However, a vari¬ 
ety of short-baseline anomalies [13-17] hint at the possible existence of short-wavelength 
oscillations, driven by one or more 0(1 eV^) mass-squared splittings that are much larger 
than the two splittings of the standard paradigm. These short-wavelength oscillations 
are purportedly already signihcant when the (Travel distance L)/(Energy E) of neutrinos 
in a beam is only ~ 1 km/GeV. Of course, they are still present at the far detector of 
any long-baseline experiment, where L/E is, say, ~ 500 km/GeV. In this work, we have 
explored the consequences of the short-wavelength oscillations - should they be real - for 
measurements at long baselines, especially measurements that probe CP violation. We hnd 
that these consequences could be considerable. For example, it is possible for long-baseline 
results, interpreted without taking the short-wavelength oscillations into account, to imply 
that CP violation is very small or totally absent, when in reality it is very large. In ad¬ 
dition, long-baseline measurements interpreted as determining the sole oscillation-relevant 
CP-violating phase in the standard paradigm could in fact be measuring something else. 

The large splittings hinted at by the short-baseline anomalies imply the existence of 
additional, largely sterile, neutrino mass eigenstates, beyond the three of the standard 
scenario (referred to as 3+0 in what follows). These additional mass eigenstates introduce 
not only additional splittings but also additional mixing angles and phases. For simplicity, 
we restrict ourselves to the scenario, referred to as 3+1, with only one additional mass 
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eigenstate. In this scenario, there are six mixing angles, and three CP-violating phases 
that can affect oscillation. Denoting the mass eigenstates of 3+0, as usual, as 
and the additional mass eigenstate as and defining the mass-squared splittings as = 
mf — mj (i, j = 1, 2, 3, 4 &: i / j), we have, according to present data, 

6m1i ~ (5m42 ~ <^^43 >> |<5^3il rsj \5ml2\»Sml^. (1.1) 


Since the probability of an oscillation driven by a splitting Smfj is proportional to 

sin^Ajj, where Aij = 1.27 x — ^ ^ ’ when L/E ~ 500 km/GeV, the short- 
wavelength oscillations driven by the large splittings involving 1^4 will be averaged to an 
L/E-independent value by the hnite energy resolution of any realistic detector. But these 
rapid oscillations are still present and can have a major impact. 

We perform our calculations for the 3+1 scenario as manifested in the proposed Deep 
Underground Neutrino Experiment (DUNE)^. While we do not explore 3+N scenarios with 
N>1, we expect that if the consequences of having one extra neutrino for long-baseline 
measurements are substantial, those of having more than one must be substantial as well, 
since the world with one extra neutrino is in a sense a special case of that with more than 
one. 

Previous work examining the effects of sterile neutrinos at long baselines includes 
several studies of neutrino factories feeding baselines of about 3000 km - 7500 km, with 
muon energies in the range 20 GeV - 50 GeV, focussing on effects at both near and far 
detectors [21-25]. More recent work [26] includes a study of effects relevant to T2K [27] 
and a combined study [28] for T2K, MINOS [29] and reactor experiments. Additionally, 
issues having some overlap with those addressed here for DUNE have been discussed in 
[30], and, very recently, in [31]. 

In Sec. 2 we examine, in the 3+1 scenario, the probability of r'e appearance at the 
far end of a long-baseline beam of neutrinos born as . We derive an analytical expression 
for valid in vacuum, and discuss how its qualitative features change when matter 
effects, which will be large in the case of DUNE, are taken into account. Sec. 3 focuses 
on the neutrino-antineutrino asymmetries that are possible at long baselines in the 3+0 
and 3+1 scenarios, and on how the possible asymmetries in these two scenarios compare. 
Sec. 4 presents long-baseline far detector event rates, obtained by performing realistic rate 
calculations for the two scenarios. Sec. 5 briefly discusses two important implications of 
the results presented here. In Sec. 6, we summarize and discuss what conclusions could 
possibly be drawn, and what ones could not be drawn, especially concerning GP violation, 
from a given set of long-baseline oscillation results, and conclude. 


^ The inputs we use, and the corresponding references, pertain to the erstwhile Long Baseline Neutrino 
Experiment (LBNE), which has undergone a new phase of internationalisation and expansion. This has led 
to a change in the name of the experiment, to DUNE. Nonetheless, it is expected that the configuration we 
assume here vis a vis fluxes, baseline and energies will remain largely intact [18-20]. 
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2 The 3 + 1 electron appearance probability in vacuum and matter 


For CPV discovery in long baseline super-beam experiments, the electron neutrino appear¬ 
ance probability is crucial. We discuss its analytic form in vacuum for the 3 + 1 scenario 
prior to discussing the matter case. While it is the latter that is relevant for DUNE in 
particular, and other long baseline (LBL) experiments at baselines of 0(1000) km in gen¬ 
eral, the form of the vacuum expression provides a useful template for the identification of 
terms the importance of which will be accentuated by the presence of matter. 

We adopt the following parameterisation^ for the PMNS matrix in the presence of a 
sterile neutrino. 


^PMNS — Cl(^34, <534)0(024, <^24)0(014)0(023)0(^13)'^ 13 ) 0 ( 012 ). (2.1) 


Here, in general, 0{6ij,6ij) is a rotation matrix in the ij sector with associated phase 6 ij. 
For example. 


0(024, <^ 24 ) — 


/l 0 0 0 \ 

0 cos 024 0 sin 024 

0 0 10 

\^0 — sin 024 0 cos 024 / 


; 0(014) = 


/ cos 014 0 0 sin 014 \ 
0 10 0 

0 0 10 

\^—sin 014 0 0 cos On J 


etc. 


Using the standard formula for a flavour transition oscillation probability, we have, for the 
3+1 case: 


P^e = 4|0^40e4p X 0.5 

- ARe(Uf,iU:^U;^Ue2) Sin2 A 21 + 2Im{U^iU:^U;2Ue2) sin2A2i 

- ARe{Uf,iU:^U;^Ue3) sin2 A 31 + 2Im{Uf,iU:^U;^Ue3) sin 2 A 3 i 

- ARe{U^2U:2U;^Ue3) sin2 A 32 + 2Im{U^2U:2U;3Ue3) sin 2 A 32 . ( 2 . 2 ) 


In arriving at the above expression, we have only assumed (based on Eq. 1.1 above) 
that sin^ A 4 j averages out to be 0.5 at long baselines, and similarly sin2A4j averages out 
to be 0, when z = 1, 2, 3. 


^This choice, which, at hrst, appears not to be the most general one which could be made, is motivated 
by the fact that in any parameterisation, the electron neutrino 3 -|-l appearance probability in vacuum 
turns out to be dependent on only two specihc linear combinations of the three independent phases. For 
simplicity, we have thus incorporated this at the outset using the fact that the first and second row elements 
{Uei and U^i) in Up'{^\^g will not have 634 or +4 in them. 
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After substituting the values of the Uai in terms of the mixing angles, we obtain : 


Pt = - sin^ 20 ^ 


lie 


'lie 


+ (o^ sin^ — - sin^ 26*13 sin^ 26 *^g) [ cos^ 6*12 sin^ A31 + sin^ 6*12 sin^ A32] 

+ cos((5i3)6a^ sin 20j^g [ cos 26*12 sin^ A21 + sin^ A31 — sin^ A32] 

+ cos((524)bosin 20 ^g[ cos 2012 cos^ 0i3 sin^ A21 — sin^ 0i3(sin^ A31 — sin^ A32)] 

+ cos(( 5 i 3 + <524)0 sin 20j)g sin 20^g [ — ^ sin^ 20i2 cos^ 0i3 sin^ A21 
+ cos 2013(cos^ Gi2 sin^ A31 + sin^ 0i2 sin^ A32)] 

— - sin(5i3)6a^ sin 20j)g [ sin 2A21 — sin 2A31 + sin 2A32] 

+ - sin(524)6asin20^g[cos^ 0i3sin2A2i + sin^ 0 i 3 (sin 2 A 31 — sin2A32)] 

+ ^ sin( 5 i 3 + 524)0 sin 20 j^g sin 20 ^g [ cos^ 0i2 sin 2A31 + sin^ 0i2 sin 2A32] 

+ (6^0^ — ^o^ sin^ 2012 sin^ 20j^g — ^ cos^ 0i3 sin^ 20i2 sin^ 20^g) sin^ A21, 


(2.3) 


where, 


sin 20j)g = sin 20i3 sin 023 ( 2.4) 

6 = cos 013 cos 023 sin 2012 (2-5) 

sin 20^g = sin 2014 sin 024 (2-6) 

O = cos 014 cos 024. (2-7) 


Prior to proceeding, we briefly discuss the allowed ranges for the 3 + 1 mixing angles 
that we have used in our calculations. These have been obtained using the results of [32], 
which takes all available data on short-baseline (SBL) oscillations and performs a global 
fit to constrain active-sterile mixing. |17e4p is constrained by Ve and Ve disappearance 
searches, and is equal to sin^ 0 i 4 . From [32], the 99% C.L. limit can be taken to be (with 
some extrapolation, as the result is given for 95%), 

|?7e4l^ G [ 0 , 0 . 1 ], 


which gives 


014 G [ 0 , 20 °]. 


Similarly, using and neutral current disappearance searches, one can constrain 6/^4 

and ? 7 t- 4 ) which are given by, 


I C%41 = cos 014 sin 024, 16+41 = cos 014 cos 024 sin 034 . 
The constraints on these elements are found to be, 

G [0,0.03], 1[/,4|" G [0,0.3]. 


~ 4 - 


These translate to, 


024 e [0,11°], 034 G [0,31°]. 


Apropos Eq. 2.3, we note that, as expected, the vacuum appearance probability is 
independent of the 3-4 mixing angle and the associated CP phase. This important simpli¬ 
fication, however, does not carry over to the matter case, as we show below. Secondly, Eq. 
2.3 contains terms proportional to the sines and cosines of a) the 3-1-1 CP phase < 524 , and 
b) the sum (Jia -|- < 524 ). These are interference terms, involving mixing angles from both 
the 3-|-0 and the 3-1-1 sector. In particular, as can be determined by inspection, the terms 
involving the sine and cosine of the sum of < 5 i 3 and <524 can be significantly large and lead 
to appreciable changes in both the amplitude of the overall probability and the extent of 
CP violation. These contributions become all the more significant once matter effects are 
large. We discuss this in more detail below, both in this section and the next. 

The matter eigenstates bring about a dependence on all mixing angles and phases. 
Specifically, unlike the vacuum case, the 3-4 mixing angle and its associated phase are no 
longer quiescent, and the 3-1-1 electron neutrino appearance probability exhibits a signifi¬ 
cant dependence on them. This is illustrated by Pig. 1, where we have used the General 
Long Baseline Experiment Simulator (GLoBES) [33, 34] to generate the plots. The left 
panel shows the variation of with energy for no CP violation (all three Dirac CP phases 
set to zero) and four different values of the 3-4 mixing angle. While the variation due to 
034 , shown in the left panel, is not very large, the effect of varying <534 (while keeping 034 
fixed, right panel) within its allowed range is quite significant. It is also striking that these 
large effects on the probability in the presence of matter are brought about by parameters 
which are completely absent in the vacuum expression (Eq. 2.3) and which (in our chosen 
parametrization) play no role in SBL situations. They illustrate the important role matter 
plays in invoking and enhancing the effects of sterile states at long baselines. In these plots, 
all other parameters are hxed as described in the caption. 

Pig. 2 emphasizes the dependencies discussed above from a slightly more general per¬ 
spective. In the right (3-|-l) panel, the significant differences between the solid and dashed 
lines of a given colour emphasize the role played by matter, while the equally significant 
differences between the blue and red dashed (solid) lines demonstrate the important role 
played by CP violating phases at long baselines in matter (vacuum) if they are non-zero. 
Turning to the left (3-|-0) panel, we note the relatively large differences between these 
curves and their counterparts in the right panel, underlining the signihcant effects of ster¬ 
ile neutrinos at the operational baseline for DUNE. 

In snmmary. Pigs 1 and 2 demonstrate that additional CP phases related to a eV^ 
sterile sector play an important role at long baselines. Their effects are heightened by the 
presence of matter. In addition, parameters related to the sterile sector which are dormant 
at short baselines and in vacuum-like conditions are no longer inert once baselines are long 
and matter effects are important. 
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Figure 1: P^e vs in earth matter for 1300 km. Averaging has been done for induced oscillations. 

In the left panel, the effect of varying 634 within its allowed range is shown with all the CP phases kept 
equal to 0. In the right panel, we show the effect of varying CP violating phase 534 when 634 = 30°, and 
the other phases are 0. For both panels, we set 614 = 20° and 624 = 10°, and the parameters related to the 
3+0 sector at the best-fit values specified in Sec. 4. 



Figure 2: (both for vacuum and matter) for 3+0 (left panel) and 3+1 (right panel) vs. energy. The 

red curves represent the CP conserving case, while the blue ones depict the case with phases set to non-zero 
fixed values (see the plot label). For the blue curve in the left panel, the sole 3+0 phase 5cp was taken 
as 30°. Normal hierarchy is taken to be the true hierarchy here, and parameters related to the 3+0 sector 
have been set at the best-fit values specified in Sec. 4. 


3 A discussion of Neutrino-Antineutrino asymmetries in matter 

The consideration of CP violation in terms of an asymmetry defined at the probability level 
provides additional insight into the conclusions which can be reliably drawn from data if we 
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do not know whether 3+0 or 3+1 is the choice nature has made. Consider the asymmetry 
defined as, 


^a/3 ^ P{Vc^ _ _ ^-P«/3 _ , 

P{^a Vj^) + P{l'a —>■ P{l^a l^ji) + P{T^a 

We begin by noting an important difference between the 3+0 and 3+1 scenarios with 
respect to the numerator APa 0 of A'^p. In vacuum, CPT invariance implies that P{vi 3 —)■ 
Ua) = P{ya Ujs), which in turn implies that APp^ = —APap, and in particular that 
APap = 0 when jS = a. Thus, when there are only three neutrino flavors, there are only 
three independent potentially non-zero CP-violating differences AP^p to be measured: 
APefj,, APp^r and APre- Now, conservation of probability implies that for any number of 
flavors, 

P{va -i- vp) = l and ^ P{va vp) = 1, 

/3 /3 

where the sums are over all /3, including j3 = a. It follows that AP^p = 0. Then, since 
APap = 0 when P = a, we conclude that in vacuum, 

AP^p = 0. (3.2) 

p^a 

When there are only three flavors, this constraint implies that AP^p + APgr = 0 and that 
APpe + AP^r = 0. Since APp^ = —AP^p, it follows that, 

APe^ = AP^r = APre. (3.3) 

That is, the three “independent” CP-violating differences are equal. In particular, if there 
are only three flavors, it is not possible for CP invariance to hold in one oscillation channel, 

such as —>■ W, and yet be violated in another channel, such as —)• . 

This situation changes when there are more than three flavors. For e.g., when there 
are four flavors, as in the 3+1 scenario, there are six independent potentially non-zero 
differences APap: APe^, AP^r, APre, APes, AP^s and APrs, where s refers to the sterile 
flavor. Now the constraint of Eq. 3.2 gives rise only to relations like 

= APf,r + AP^s- (3.4) 

It is now perfectly possible for AP^e{= —APe^), the CP-violating difference that will be 
the first to be probed experimentally, to be zero, while the differences AP^r and AP^g in 
other oscillation channels that are challenging to study, are large^. 

In Fig. 3, we show the spread of at L = 1300 km for cases chosen to illustrate 
some of the important features that arise due to the presence of a fourth, sterile state. The 

®We note that any long baseline experiment involves earth-matter effects, which break CPT (in addition 
to CP). Such breaking is extrinsic, and due to the asymmetry of the earth matter through which the 
neutrinos propagate. While this may appear to destroy the conclusions reached above, which depend on 
CPT invariance, this is not the case as long as an experiment seeks to measure intrinsic (i.e.driven by 
phases in the mixing matrix) CP violation and devises appropriate means to do so. 

^Henceforth we drop the superscripts a and /3 and take to denote the asymmetry for a = fj, and 
P = e. 
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left-hand panels were created with all CP-violating phases set to zero, so the asymmetries 
shown in these panels are from matter effects only. The right-hand panels were created 
allowing the sole 3-1-0 CP phase 5cp to vary over its entire physical range in the case of 
3-1-0, and the three CP phases < 5 i 3 , ^24 and ^34 to vary over their entire ranges in the case 
of 3-|-l. Thus, these panels show the impact of intrinsic CP violation. In all panels, the 
red curve(s) are for the 3-1-0 case, and the blue ones for the 3-1-1 case. The top two panels 
assume a normal hierarchy, and the bottom two an inverted hierarchy. In creating all 
panels, the mass splittings and mixing angles of the 3-1-0 sector were set to the best-fit 
values specified in Sec. 4, the splitting Am|^ of the 3-1-1 sector was set to 1 eV^, and the 
3-|-l mixing angles On, O 24 and ^34 were varied over their allowed ranges. 

With one exception (see below), to create the curves in each panel of Fig. 3 for each of 
the two scenarios, 3-1-0 and 3-1-1, we varied the corresponding parameters until we found the 
parameter set that maximizes (minimizes) the energy-integrated asymmetry Ai,p for that 
scenario. The energy-dependent asymmetry was then plotted vs. energy for this parameter 
set as a solid (dashed) curve. (Note that since it is the energy-integrated asymmetry that 
is being extremized, it is possible for the 3-1-0 energy-dependent asymmetry to be more 
extreme than that for 3-1-1 for a limited range of energy, despite the fact that the 3-1-0 
scenario is in a sense, a special case of 3-1-1.) The one exception to our procedure is that, 
since the 3-1-0 sector parameters other than (5cp were held fixed throughout, in creating the 
left-hand panels, no 3-1-0 parameters were varied, so there is only a single curve, shown as 
solid, for 3-1-0. 

From the left-hand panels of Fig. 3, we see that when CP is conserved, the neutrino- 
antineutrino asymmetry vs. energy is quite similar in the 3-1-0 and 3-1-1 scenarios. In the 
3-|-l scenario, this asymmetry is confined to a rather narrow band as the 3-1-1 mixing angles 
are varied. Although it is not shown, we have found that it is confined to a similar narrow 
band in the 3-1-0 scenario as the 3-1-0 parameters 023 and are varied within their 

experimental uncertainties. Clearly, if an experiment were to measure an asymmetry vs. 
energy that consistently lies outside the similar, narrow 3-1-0 and 3-1-1 bands that correspond 
to CP conservation, we would have evidence that CP is violated so long as nature has chosen 
either the 3-1-0 or 3-1-1 scenario. However, a measured asymmetry between —)• and 
that lies within these similar narrow bands would not unambiguously signal that 
there is no CP violation in neutrino oscillation. As explained above, when there are more 
than three flavors, as in the 3-1-1 scenario, it is possible for there to be little or no CP 
violation in one oscillation channel, and yet a large CP violation in some other channel. In 
addition, for either the 3-1-0 or 3-1-1 case, it might happen that for some non-zero values of 
the CP-violating phases and mixing angles slightly different from those corresponding to 
the CP-conserving bands, the asymmetry still lies within those bands within uncertainties. 

The right-hand panels in Fig. 3 show that when intrinsic CP is violated, Ai,p can be 
anywhere in a large range. Moreover, for 3-1-1, this range is much larger than for 3-1-0, 
and includes almost all of the 3-1-0 range. Thus, we see that sterile neutrinos with 0(1) 
eV^ masses can very substantially impact CP-violation measurements at long baselines. 
While a measured asymmetry outside the band allowed for 3-1-0 would be evidence for new 
physics beyond 3-1-0, one inside that band would leave uncertain the precise origin of the 
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observed CP violation. 

To understand why the effect of a fourth, sterile neutrino on CP violation at long base¬ 
lines can be so large, one notes that CP-violating phases affect physics through interferences 
between amplitudes. As pointed out in [26], around the first maximum of the atmospheric- 
wavelength oscillation, where the long-baseline experiments work, or will work, the (new, 
short wavelength oseillation) - (atmospherie-wavelength oscillation) interference, and the 
(atmospheric-wavelength oscillation) - (solar-wavelength oscillation) interference, can eas¬ 
ily be of comparable size. Then, if the CP phases are right, 3-1-1 can be quite different 
from 3-1-0. 

In the next section, we see how the probability-level results of this section and the 
previous one translate into observable consequences for DUNE by calculating event rates. 



Figure 3 : The neutrino-antineutrino asymmetry A^v vs. energy E. See text for explanation and discus¬ 
sion. 

4 Event Rates at DUNE in the 3+1 and 3+0 scenarios 

Having discussed some of the salient features of the probability both in vacuum and 
matter at baselines characteristic of DUNE, we now perform event rate calculations by 
realistically simulating the experiment. 

We recall that DUNE (with specifications very similar to LBNE,[18], [19]) is a proposed 
future super-beam experiment, to be located in the United States, with a main aim of 
establishing or refuting the existence of CPV in the leptonic sector. In addition to this 
primary goal, this facility will also be able to resolve other important issues like the mass 
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hierarchy and the octant of 023- The super-beam will originate at the Fermilab. The 

primary beam simulation assumes a 1.2 MW - 120 GeV proton beam which will deliver 
10^^ protons-on-target (POT) per year. A 35-40 kt Liquid Argon (LAr) far-detector will 
be housed at the Sanford Underground Research Facility in the Homestake mine in South 
Dakota, 1300 km away. The experiment plans to have a total of 10 years of running, 
divided equally between neutrinos and anti-neutrinos, corresponding to a total exposure 
of 35 X 10^^ kt-POT-yr. The other experimental details such as signal and background 
definitions as well as the detector efficiencies taken in this work are the same as those in 
[20], except with the difference that we have not considered tau events in the backgrounds. 
The detector efficiencies for both and P^g events are close to 80% with somewhat less 
efficiency for P^g. 

In order to facilitate the drawing of physics conclusions, we assume certain values and 
ranges for neutrino oscillation parameters in the standard 3-flavour paradigm, which are 
motivated by their current measured ranges and best fit values. Specifically, 


• 012 and 013 have been fixed at 33.48° and 8.5° respectively [35]. 


• We assume that the 2-3 mixing is near-maximal i.e. 023 = 45°.^ 


• We fix at 7.5 x 10“® eV^ and at 2.457 x 10“^ eV^ (—2.374 x 10“^ eV^) 

for NH (IH). These particular choices have been taken from the analysis of global 
data [35]. 


We note that as stated above, standard global analyses ([36], [37], [35]) assume the 3-1-0 
scenario, but, as has been demonstrated in e.g. [38], their conclusions remain very robust 
in the presence of sterile neutrinos. Our assumed ranges for the sterile sector mixing angles 
corresponding to the 3-|-l scenario are as stated in Sec. 2. In addition, we assume Am|^ 
to be 1 eV^, and vary dis, <524 and <534 for 3-I-I and dcp for 3-1-0 in the entire possible range 
of [—180°, 180°]. Finally, we use the fluxes provided in [39]. 


® We note that the currently-allowed Scr range on 023 is [38.3°, 53.3°] with the best fit at 42.3°(49.5)° 
for NH (IH) [35]. The 023 best fit values from the global analyses [36, 37] are different from [35]. However, 
in this work, we make the simplifying assumption of maximal mixing, relegating a more rigorous statistical 
analysis to a follow-up work. 
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Figure 4: Neutrino and anti-neutrino event rates in DUNE plotted as a function of the reconstructed 
neutrino energy. The vertical spread for a given color for an energy bin shows the maximum and the 
minimum events rates possible. 


Fig. 4 shows the spread of binned events as a function of the reconstructed neutrino 
energy for the 3+0 case and the 3+1 case. For the 3+0 case, we varied only (icp in the 
range [—180°, 180°] to obtain the events band shown in red. For 3+1, we chose three sets 
of 014, 024 values - (20°, 10°), (15°, 10°) and (5°, 5°). For all these three cases, we varied 034 
in the range [0,30°] and the phases ( 5 i 3 , 624 and ^34 in the range [—180°, 180°] each. The 
resulting event-bands are shown in blue, green and magenta, respectively. The left (right) 
panels show the neutrino (anti-neutrino) rates, while the top (bottom) panels are for the 
NH (IH) scenario. 

It can be seen that for all three sets of 0 i 4 , 024, the 3+1 band can potentially encompass 
the 3+0 band, leading to substantial degeneracy. When the number of events falls in the 
overlapping region between these two bands (which is the red region in Fig. 4), there is 
considerable ambiguity as to whether the events are produced by a certain value of (5cp in 
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the 3+0 sector or by some combination of @ 34 , ( 5 i 3 , 624 and ^34 in the 3+1 sector. 

Fig. 4 also shows that the 3+1 band gets wider as the values of 014,024 and hence the 
effective mixing angle sin20^g increase. Indeed, for sufficiently large 3+1 mixing angles, 
the 3+1 band is substantially larger than its 3+0 counterpart. An observed surfeit or a 
dearth of events compared to those expected in the 3+0 case, especially near the event 
maxima (around the region 2-4 GeV), could be a pointer to the presence of sterile states. 


5 Implications 

Prior to summarizing our conclusions, it is useful to discuss certain implications which 
arise from our results. We have noted above that in the presence of even a single sterile 
neutrino, conclusions such as, a) Whether CP is conserved or violated, and b) if the latter, 
whether the violation is ascribable to the active neutrinos or the additional sterile neutrino, 
or a combination of the two, are all rendered significantly ambiguous. An important 
consequence is thus the need for an improved synergistic linkage between the global LBL 
and SBL efforts [40], since it appears that results obtained in the former cannot be correctly 
interpreted without definitive conclusions drawn from the latter. 

Our work also has ramifications for near detector (ND) design and physics requirements 
in DUNE and other LBL experiments. Even if there is no sterile sector and the standard 
three-family scenario is nature’s choice, in order to fully exploit the CPV capabilities of a 
far detector (ED), the ND must establish the expected number of events at the ED in the 
absence of oscillations with very high precision, in order to ensure systematic errors stay well 
below statistical ones. The rates expected at the ED depend on fluxes and cross-sections 
measured, along with their energy dependence, to significantly high accuracy at the ND for 
all four species of neutrinos, I'e, In the 3 + 0 scenario, these measurements, while 

very demanding, are assumed to be made under conditions where there are no oscillations 
between the source and the ND. This task is rendered significantly more complex, however, 
in the presence of a sterile sector capable of altering the fluxes between the source and the 
ND over the planned distance of ~ 500 m in DUNE. Historically, uncertainties in source 
fluxes and cross-sections have always been a major limiting factor for neutrino experiments 
seeking precision in oscillation studies. The high intensity of the DUNE beam and the 
presence of a highly capable and precise ND, designed to overcome these limitations, are 
thus rendered even more crucial than assumed earlier in order to deconvolute the added 
complexities arising due to a sterile sector, if it were to be present in nature. Moreover, we 
note that current data allow the mass-squared splitting in the sterile sector to be sufficiently 
large such that oscillations due to it can average out even at the ND distance of ~500 m. 
Such an eventuality further underscores the need for a strong connection between the SBL 
and LBL experimental programs, since these oscillations, even though not directly observed 
at the ND, would nonetheless affect its crucial service task of accurate flux determination. 
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6 Summary and Conclusions 


To summarize, we have studied the effects of the additional mixing angles and CP phases 
in the case of a 3+1 sterile sector on the determination and measurement of CPV at 
long baselines for the DUNE experiment. From a probability analysis, we show that the 
effects of the additional CP phases can be large at its chosen baseline of 1300 km. These 
effects, which arise from large interference terms (between the 3+0 and 3+1 sectors) in 
the appearance probability, are accentuated by the presence of matter, which additionally 
brings in contributions from sterile-sector mixings and phases which are dormant at short 
baselines. From event rate calculations, we show that the presence of a sterile sector 
manifests itself in measurably altered rates in energy bins across the spectrum, without 
signihcant distortion in the shape. This alteration in event rates increases, as expected, for 
larger values of the mixing angles connecting the active and sterile sectors. 

Importantly, the presence of a sterile sector obfuscates conclusive determinations of 
CP violation or conservation at the far detector, and makes uncertain the ability to ascribe 
any perceived CPV to a unique phase in the 3+0 sector. Thus, the linkage between the 
presently planned long and short baseline programs must be explored and strengthened. 
Until the presence of an ~ eV^ sector is conclusively ruled out, our work emphasizes 
the need for a complementary SBL sterile-search program and for a highly capable and 
versatile near detector for DUNE, enabling it to reduce systematics to low levels so that it 
may achieve its stated primary goals for CPV detection. 

Acknowledgements 

We thank Mary Bishai for help with DUNE fluxes and for patiently answering questions 
related to DUNE. RG and BK thank William Louis for very useful discussions. RG ac¬ 
knowledges useful discussions with Sandhya Ghoubey, Amol Dighe and S. Uma Sankar. 
MM thanks Animesh Chatterjee for very helpful discussions. RG and BK are grateful to 
the Kavli Institute for Theoretical Physics for support and hospitality. RG acknowledges 
support from the Fermilab Neutrino Division and the DUNE Project at Fermilab while 
this work was in progress in the form of an Intensity Frontier Fellowship in 2014 and a 
Guest Scientist position in 2015. RG, MM and SP also acknowledge support from the XII 
Plan Neutrino Project of the Indian DAE. BK thanks the Mainz Institute for Theoreti¬ 
cal Physics for its hospitality and support while this paper was completed. This research 
was supported in part by the US National Science Foundation under Grant No. NSF 
PHYl 1-25915. Fermilab is operated by Fermi Research Alliance, LLC under Gontract No. 
DE-AC02-07CH11359 with the US Department of Energy. 

References 

[1] M. Thomson, talk at WINP at Brookhaven, February 2015. 

’’https://indico.bnl.gov/conferenceDisplay.py?confId=918” (2015). 

[2] A. Norman, AIP Conference Proceedings 1666, 110001 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915573. 


- 13 - 


[3] C. Farnese, AIP Conference Proceedings 1666, 110002 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915574. 

[4] S. Dusini, AIP Conference Proceedings 1666, 110003 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915575. 

[5] A. B. Sousa, AIP Conference Proceedings 1666, 110004 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915576. 

[6] A. K. Ichikawa, AIP Conference Proceedings 1666, 130001 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915579. 

[7] F. J. P. Soler, AIP Conference Proceedings 1666, 130002 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915580. 

[8] Y. Hayato, AIP Conference Proceedings 1666, 130003 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915581. 

[9] R. J. Wilson, AIP Conference Proceedings 1666, 130004 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915582. 

[10] T. Patzak, AIP Conference Proceedings 1666, 130005 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915583. 

[11] K. Long, AIP Conference Proceedings 1666, 130006 (2015), URL 

http://scitation.aip.org/content/aip/proceeding/aipcp/10.1063/1.4915584. 

[12] A. de Gouvea et al. (Intensity Frontier Neutrino Working Group) (2013), 1310.4340. 

[13] A. Aguilar-Arevalo et al. (LSND), Phys.Rev. D64, 112007 (2001), hep-ex/0104049. 

[14] A. Aguilar-Arevalo et al. (MiniBooNE), Phys.Rev.Lett. 102, 101802 (2009), 0812.2243. 

[15] G. Mention, M. Fechner, T. Lasserre, T. Mueller, D. Lhuillier, et al., Phys.Rev. D83, 073006 
(2011), 1101.2755. 

[16] T. Mueller, D. Lhuillier, M. Fallot, A. Letourneau, S. Cormon, et al., Phys.Rev. C83, 054615 
(2011), 1101.2663. 

[17] A. Aguilar-Arevalo et al. (MiniBooNE), Phys.Rev.Lett. 110, 161801 (2013), 1207.4809. 

[18] J. Hewett, H. Weerts, R. Brock, J. Butler, B. Casey, et al. (2012), 1205.2671. 

[19] C. Adams et al. (LBNE Collaboration) (2013), 1307.7335. 

[20] M. Bass et al. (LBNE), Phys. Rev. D91, 052015 (2015), 1311.0212. 

[21] A. Donini, M. Lusignoli, and D. Meloni, Nucl.Phys. B624, 405 (2002), hep-ph/0107231. 

[22] A. Dighe and S. Ray, Phys.Rev. D76, 113001 (2007), 0709.0383. 

[23] A. Donini, K.-i. Fuki, J. Lopez-Pavon, D. Meloni, and O. Yasuda, JHEP 0908, 041 (2009), 
0812.3703. 

[24] O. Yasuda, pp. 300-313 (2010), 1004.2388. 

[25] D. Meloni, J. Tang, and W. Winter, Phys.Rev. D82, 093008 (2010), 1007.2419. 

[26] N. Klop and A. Palazzo, Phys. Rev. D91, 073017 (2015), 1412.7524. 

[27] K. Abe et al. (T2K), Phys.Rev.Lett. 112, 061802 (2014), 1311.4750. 

[28] B. Bhattacharya, A. M. Thalapillil, and C. E. Wagner, Phys.Rev. D85, 073004 (2012), 
1111.4225. 


- 14 - 


[29] P. Adamson et al. (MINOS), Phys.Rev.Lett. 112, 191801 (2014), 1403.0867. 

[30] D. Hollander and I. Mocioiu, Phys.Rev. D91, 013002 (2015), 1408.1749. 

[31] J. M. Berryman, A. de Gouvea, K. J. Kelly, and A. Kobach (2015), 1507.03986. 

[32] J. Kopp, P. A. N. Machado, M. Maltoni, and T. Schwetz, JHEP 1305, 050 (2013), 
1303.3011. 

[33] P. Huber, M. Lindner, and W. Winter, Comput.Phys.Commun. 167, 195 (2005), 
hep-ph/0407333. 

[34] P. Huber, J. Kopp, M. Lindner, M. Rolinec, and W. Winter, Comput. Phys. Commun. 177, 
432 (2007), hep-ph/0701187. 

[35] M. Gonzalez-Garcia, M. Maltoni, and T. Schwetz, JHEP 1411, 052 (2014), 1409.5439. 

[36] F. Capozzi, G. Fogli, E. Lisi, A. Marrone, D. Montanino, et al., Phys.Rev. D89, 093018 
(2014), 1312.2878. 

[37] D. Forero, M. Tortola, and J. Valle, Phys.Rev. D90, 093006 (2014), 1405.7540. 

[38] A. Esmaili, E. Kemp, O. Peres, and Z. Tabriz!, Phys.Rev. D88, 073012 (2013), 1308.6218. 

[39] Mary Bishai, private communication (2014). 

[40] Joseph Zennamo, 48th Annual Fermilab User’s Meeting (June 2015), URL 

https://indico.fnal.gov/getFile.py/access?contribId=54&sessionId=23&resId= 
0&materialId=slides&confId=8982. 


- 15 - 


